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摘 要 : 为 探 明 各 气象 因素 对 锡林郭勒 草原 的 参考 作物 蒸 散 量 影响 及 对 参考 作物 蒸 散 量 预 


勒 草 原 的 6 个 国家 级 地 了 


排序 分 别 建 立 基 于 多 变量 时 间 序 列 CAR(Controlled Auto-regressive) 的 CAR-ETo 模 型 。 结 果 表 明 : (1 


测 ,利用 内 蒙古 锡 林 郭 


气象 站 点 获取 的 气象 数据 与 其 计算 得 到 的 PM-ET, 进 行 关联 度 分 析 ,按照 关联 系数 从 大 小 


各 气象 因子 与 


PM-ET 关 联系 数 呈 现 由 荒漠 草原 区 -典型 草原 区 - 草 旬 草原 区 过 滤 减 小 的 分 布 ,日 最 低 气 温 .平均 气温 与 PM-ET 关 


联 度 最 小 ,平均 日 照 时 数 与 PM-ET 关 联 度 最 大 ,为 0.7293。(2) 处 于 草 甸 草 原 区 .荒漠 草原 区 


两 因子 建立 的 CAR-ET, 精 度 较 高 
型 预测 精度 总 体高 于 HS-ET, 与 P 


;典型 草原 区 还 需 风速 为 模型 第 三 因子 ,(3) 通过 模型 预测 


.通过 输入 日 照 与 气温 
精度 验证 ,CAR-ET。 模 


MT-ET 方 法 ,结合 关联 度 分 析 可 得 出 对 当地 参考 作物 蒸 散 和 


是 有 显著 影响 的 气象 因 


子 , 可 为 锡林郭勒 草原 气象 监测 布设 方案 .参考 作物 蒸 散 量 确定 与 草原 生态 恢复 提供 理论 依据 。 
关键 词 : SS WAR; 灰色 关联 度 分 析 ; 多 变量 时 间 序列 模型 ; 锡林郭勒 草原 


参考 作物 蒸 散 量 (ETu) 是 计算 作物 需 水 量 的 重 
要 参数 ,准确 计算 ET 对 区 域 水 资源 评价 与 高 效 利 
用 意义 重大 ,ET 在 水 文 模拟 灌溉 制度 的 制定 水 
量 平衡 计算 .气候 干 湿 等 分 析 中 有 着 重要 的 作用 ， 
是 农田 水 分 管理 与 灌溉 决策 的 重要 参数 ,也 是 评价 
气候 干 湿 程 度 、 植 被 耗 水 量 的 重要 指标 之 一 。 目 
前 ,常用 的 ET 计算 方法 为 FAO-56 Penman-Monteith 
(FA056-PM) 法 中 简称 (PM-ET。) 法 。PM-ET, 模 型 需 
要 太阳 净 辐 射 (R,) 、 最 高 和 最 低 气 温 (7,、7,,) 、 相 
对 湿度 (RA) 和 2.0 m 高 度 风速 (w) 等 数据 。 但 在 类 
似 锡林郭勒 草原 这 样 的 地 域 广 阔 、 能 够 提供 E76 计 
算 所 需要 的 全 部 气象 因子 的 国家 气象 站 点 较 少 ,就 
要 求 输入 气象 因子 较 少 的 ET 估算 方法 。 例 如 Har- 
greaves-Samani 方程 “简称 HS-ETu, 仅 需 最 高 和 最 
AR AT Yi (Tas Ts) 就 可 估算 E760 ,同时 Penman- 
Monteith 温度 (PMT) 法 也 仅 需 气温 数据 ” ,以 上 2 种 
基于 气温 数据 估算 参 考 作物 蒸 散 量 的 方法 在 研究 
中 得 到 广泛 应 用 ””。 近 年 来 , 随 着 计算 机 技术 不 断 
发 展 产 生 了 许多 新 的 估算 方法 ,如 Ferreira 等 "使 用 


收 稿 日 期 : 2021-03-01; 修订 日 期 : 2021-09-01 


机 器 学 习 与 深度 学 习 算 法 利用 小 时 数据 估算 参考 
作物 蒸 散 量 , Yamac 等 "利用 3 种 不 同 的 机 器 学 习 
算法 和 4 种 气象 资料 情景 估算 马 铃 蔚 作物 的 日 燕 散 
He, 4E T RSE |_A ey HE Copulla KAU XT Mn] ET, 
的 气象 因素 进行 联合 分 布 构建 ,揭示 不 同 变量 间 的 
相关 结构 ,建立 多 元 气象 因素 对 Tv 的 联合 分 布 模 
型 。 虽 然 这 类 模型 不 需要 复杂 物理 参数 ,但 是 这 些 
预测 模型 对 建 模 时 间 序 列 数据 有 平稳 性 .独立 性 与 
正 态 性 等 要 求 。 多 变量 时 间 序 列 CAR (Controlled 
Auto-Regressive) 模 型 结合 了 一 维 时 间 序 列 分 析 和 
回归 分 析 2 类 数理 统计 方法 的 优点 ,不 仅 考虑 到 事 
物 发 展 的 自身 运动 规律 ,顾及 了 环境 因子 的 作用 ， 
可 以 有 效 避 免 处 理 复杂 时 间 序 列 问题 而 引起 的 误 
差 ” ,多 变量 时 间 序 列 模型 如 今 在 水 文 预报 地 下 水 
埋 深 呈 ,气象 “等 得 到 广泛 应 用 , 且 效 果 显 著 。 

本 文通 过 将 PM-ET 模 型 计算 值 作为 标准 值 , 利 
用 锡林郭勒 草原 地 面 气象 站 所 得 气象 因子 结合 ; 
联 度 分 析 建 立 多 变量 时 间 序 列 CAR-ET 模 型 并 验证 
预测 精度 ,与 其 他 常见 的 ET 计算 方法 进行 比较 , 提 
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出 适合 锡林郭勒 草原 的 7, 计 算 方法 ,以 期 为 锡 林 
郭 勒 草原 设置 气象 监测 .计算 参考 作物 蒸 散 量 与 草 
原生 态 恢 复 提供 理论 指导 。 


1 数据 与 方法 


1.1 研究 区 概况 

锡林郭勒 位 于 内 蒙古 中 东部 ,蒙古 高 原 南 部 。 
东西 长 约 700 km, 南 北 宽 500 km, 总 面积 约 为 2.03x 
10 km ,其 中 天 然 草场 面积 占 总 面积 95% , 约 为 
内 蒙古 草原 面积 的 1/4 ,是 中 国 重 要 的 畜产 品 基地 与 
生态 保护 屏障 号 。 地 处 中 纬度 西风 带 , 属 中 温带 半 
干旱 .干旱 大 陆 性 季风 气候 ,其 具有 风 大 ,干旱 , 寒 
冷 的 气候 特点 ,土壤 以 粉 砂 土 为 主 下 ,地 势 南 高 北 
低 , 东 、 南 部 多 低 山 夺 陵 , 西 .北部 地 形 平 坦 。 锡 林 
郭 勒 根据 牧区 自然 地 理 条 件 差异 可 分 为 3 个 地 带 性 
植被 类 型 区 域 ,东部 以 草 旬 草原 区 为 主 ,中 部 以 典 
型 草原 区 为 主 ,西部 以 荒漠 草原 区 为 主 ”1。 
1.2 数据 来 源 

锡林郭勒 盟 气 象 局 共 设 15 个 旗 ( 县 市、 区) 级 
气象 局 ( 台 、 站 ), 但 只 有 国家 级 地 面 站 提供 的 数据 
才能 满足 PM-ET 方 法 计算 参考 作物 蒸 散 量 。 本 文 
所 用 到 的 气象 数据 是 在 国家 气象 科学 数据 中 心 获 
取 , 包 括 了 锡林郭勒 盟 6 个 国家 级 地 面 站 。 所 用 和气 
象 因子 数据 为 日 最 高 气温 .日 最 低 气 温 .平均 相对 
湿度 日照 时 数 .平均 风速 ,研究 区 气象 ( 台 、 站 ) 点 
分 布 如 图 1 所 示 , 各 站 点 1980 一 2019 年 平均 气候 资 
料 汇总 如 表 1 所 示 。 
1.3 参考 作物 蒸 散 量 计算 模型 
1.3.1 FAO-56 Penman-Monteith(PM-ET,) 法 其 表 
达 式 为 : 


900 
0.408A(R, -G)+y u,(e, —e,) 
= T +273 ( 1 ) 
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图 1 研究 区 站 点 分 布 


Fig. 1 Meteorological station distribution map 


式 中 :7 为 参考 作物 蒸 散 量 ,mmd';A 为 温度 随 饮 
和 水 气压 变化 的 斜率 ,kPa…°C7';RR, 为 作物 表面 的 净 
辐射 , MJ.m*…d ;6G 为 土壤 热 通 量 , MJ md ;7 为 
干 湿 计 常数 ,kPa*C7 ;7 为 日 平均 气温 ,"C; 岂 为 2m 
高 度 处 风速 ,ms'!;e, 为 饱和 水 气压 ,kPa;e, 为 实际 
水 汽 压 kPa。 
1.3.2 多 变量 时 间 序 列 (CAR) 模 型 ， 多 变量 自 回归 
模型 建 模 方法 主要 有 2 种 :一 是 采用 一 般 最 小 二 乘 
法 建立 回归 模型 的 方法 ;二 是 递 推 最 小 二 乘法 进行 
模型 参数 估计 的 方法 。 在 DPS 数 据 人 处 理 系 统 ” 中 
采用 的 是 第 2 种 方法 建立 CAR 模 型 。 

为 了 不 失 一 般 性 ,假设 用 mm 个 变量 的 时 间 序 列 
组 建 n 阶 CAR 模 型 ,其 形式 为 : 

Y=ay FV 2 ya FO 9%, +O H+ 


batia Din Xin + Big Xp, T bakat 
Uh 2s" ` DMs ia 本 ss EO 0 + 


(2) 


bm i tba. "Ob Nate, 
对 模型 采用 递 推 最 小 二 乘法 进行 参数 估计 , 令 
0=[a,, ao 
Days bais batts Dy] 
z E[n Jen Mia Mics Migs Bias 


ee Mons Xo 1-09 Xo p19 2 -29 zery: Airal 


(3) 


T 


as Ay +0341, 则 CAR 模 型 可 以 表示 为 ; 
#1 1980 一 2019 年 各 站 点 气候 资料 
Tab.1 The average climatic data of each station from 1980 to 2019 were summarized 
气象 站 点 生态 地 理 分 区 ”日 最 高 气温 /CC ”日 最 低 气 温 /C ”平均 相对 湿度 /% 日 照 时 数 /h 平均 风速 /ms 

阿 巴 嘎 旗 典型 草原 区 9.252 -4.001 55.072 8.142 3.166 
RARI ESR 9.257 -4.213 60.876 8.157 2.852 
西 乌 珠 称 草 多 草原 区 8.928 -4.194 58.485 7.681 3.334 
多 伦 典型 草原 区 10.041 -3.308 59.663 8.038 3.172 
二 连 浩特 RJRK 12.229 -1.623 46.672 8.715 3.850 
锡林浩特 典型 草原 区 10.295 -3.036 55.503 8.132 3.343 
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y=z 0+e, (4) 
式 中 :7 表示 和 矩阵 的 转 置 。 根 据 该 式 ,在 时 刻 /,6 的 
递 推 最 小 二 乘 估计 值 为 6) : 
6=0 +KO -2z70 1) 
Piz 
a P E (5) 
P,= t0 -KP,, 
可 利用 和 N 组 观察 值 所 得 到 的 CAR(n) 模 型 来 计 
算 残 差 平 方 和 5S(n): 
Sn)= > e (6) 


Kt= 


NPR: 

ê =y, -z ô, (7) 
式 中 :ao 为 遗忘 因子 ,一 般 取 0.9~1.0, 如 果 要 求 遗忘 
历史 数据 的 速度 快 些 ,那么 应 趋 近 于 0.9; 如 果 要 
求 遗忘 历史 数据 的 速度 慢 些 ,那么 应 趋 近 于 1.0; 
如 果 a = |, 那么 历史 数据 永远 不 会 遗忘 掉 。 

XT CAR 模型 的 定 阶 是 根据 已 知 的 V 个 样本 (x， 
yist=1,2,°°°,Nsj=1,2,°++,m), ARB ES ee 
地 对 系统 拟 合 CAR 模型 ,并 且 依 次 对 相 邻 的 2 个 
CAR 模 型 采用 下 检验 的 方法 判断 模型 的 阶 次 增加 是 
否 合适 。 

对 CAR(o) 模 型 中 的 某 些 系数 是 否 为 0 进行 下 检 

验 , 以 决定 模型 的 真实 阶 及 其 时 清 ,从 而 得 到 真实 
模型 参数 估计。 
1.3.3 灰色 关联 度 分 析 灰色 关联 度 分 析 与 数理 统 
计 学 的 相关 分 析 不 同 ,关联 度 分 析 是 基于 灰色 系统 
的 灰色 过 程 ,分析 是 因素 间 时 间 序 列 的 比较 ,同时 
关联 分 析 不 要 求 数据 太 多 ,主要 研究 动态 过 程 ,而 
相关 分 析 则 以 静态 研究 为 主 ,因此 ,关联 度 分 析 适 
应 性 更 广 。 关 联 度 分 析 在 农业 ,水文 .气象 等 领域 
十 分 广泛 ,如 路 璐 等 3 运用 关联 分 析 法 分 析 了 水 文 
气象 因子 对 土壤 水 分 变化 起 到 主要 影响 与 次 要 影 
响 , 王 秋 京 等 中 利用 灰色 关联 分 析 方 法 分 析 了 黑龙 
江 省 农业 气象 灾害 演变 特征 , 刘 锋 等 提出 了 改进 
灰色 关联 从 整体 角度 寻求 系统 的 最 优 设计 方案 。 
1.3.4 模型 评价 指标 ”通过 CAR-ETo 模 型 与 PM-ET,。 
方法 计算 出 的 结果 比较 ,采用 决定 系数 ( 尼 ) E 
方 根 误差 (RMSE) Nash- Sutcliffe 建 模 效率 系数 
(NSE) 对 模型 模拟 效果 进行 评价 : 

公式 如 下 : 


R= (8) 
$o -0 Se.-a] 
Seo 
RMSE =| =— (9) 
QL -CF 
NSE = | - =) (10) 
XQ..- 2.) 


式 中 : Q! Q, 分别 为 预测 值 与 观察 值 (4=1,2,3,*…， 
n) 之 间 的 差异 ; 0, SQ, 表示 预测 值 与 观察 值 的 
均值 。 
1.4 数据 分 析 

利用 Excel 2016 软件 按照 每 1 行为 1 个 样本 ， 
1 列 为 1 个 变量 整理 数据 ,将 输入 因子 放置 在 左边 若 
干 列 ,将 右边 最 后 一 列 设 为 输出 因子 ,DPS7.05(Da- 
ta processing system) 数 据 处 理 将 竺 分 析 数 据 按照 规 
定格 式 定 义 为 数据 块 , 分 别 在 灰色 系统 方法 与 时 间 
序列 下 模块 下 进行 灰色 关联 度 分 析 与 建立 多 变量 
时 间 序 列 模型 ,应 用 QGIS 3.14 软件 的 反 距 离 加 权 法 
插值 将 各 个 气象 站 点 得 到 的 灰色 关联 系数 插值 在 
研究 区 ,并 利用 Origin 2018 软件 作 图 。 


2 结果 分 析 


21 参考 作物 蒸 散 量 与 各 气象 因子 关联 度 分 析 
根据 关联 度 分 析 得 到 的 多 年 平均 PM-ET, 在 各 
站 点 与 各 参数 的 关联 系数 图 如 图 2 所 示 。 多 年 平均 
PM-ET, 与 各 气象 因子 的 相关 联 程度 受 参 数 本 映 及 
地 理 位 置 因素 的 影响 ,关联 系数 分 布 的 总 体 趋势 为 
位 于 研究 区 西北 方向 的 关联 系数 较 东 南方 向 的 关 
联系 数 更 大 ,在 研究 区 呈现 由 西 到 东 减 小 的 分 布 。 
研究 区 内 与 多 年 平均 PM-ET, 关 联 度 最 小 的 2 个 因 
子 为 多 年 最 低 气 温 的 平均 值 多 年 平均 气温 ,最 大 
为 多 年 平均 日 照 时 数 。 由 于 处 在 西部 的 二 连 浩特 
市 受 蒙古 高 压 影响 ,属于 中 温带 大 陆 性 季风 气候 、 
干旱 荒漠 草原 气候 , 蒸 散 发 潜力 强 , 降 水量 少 ” , 参 
考 作物 蒸 散 量 易 受 太阳 辐射 影响 ,处 在 研究 区 东北 
方向 的 东 马 珠 穆 沁 旗 受 东南 海洋 暧 湿 气 流 与 西北 
干 寒 气流 影响 , 属 寒 温带 大 陆 性 季风 气候 ,大 陆 型 


202112.00102v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


6 汉 壮 壮 等 :基于 多 变量 时 间 序 列 模 型 的 锡林郭勒 草原 参考 作物 蒸 散 量 预 测 1653 
(a) 多 年 平均 气温 (b) 多 年 平均 最 高 气温 (c) 多 年 平均 最 低 气温 
(d) 多 年 平均 风速 (e) 多 年 平均 空气 湿度 O 多 年 平均 日 照 时 数 


we a” 


Œ 0.424~0.496 © 0.712~0.784 


关联 系数 
国 <0.352 


>z 


= 0.496~0.568 E 0.784~0.856 
0.568~0.64 
= 0.352~0.424 E 0.64~0.712 


E 0.856~0.928 0 175 km 
E >0.928 上 


图 2 各 站 点 的 PM-ET, 与 各 参数 的 关联 系数 分 布 


Fig. 2 The distribution diagram of PM-ET, and correlation coefficient of each parameter at each station 


高 山 气 候 区 , AES BEER?! ,日 照 时 间 相 对 较 小 ， 
其 参考 作物 蒸 散 量 与 平均 日 照 时 数 关联 度 也 较 大 ， 
多 年 平均 日 照 时 数 与 PM-ET, 的 关联 系数 最 大 值 出 
现在 二 连 浩特 与 东 乌 珠 穆 沁 旗 , 关 联系 数 分 别 为 
0.6973 .0.7293 。 

2.2 不 同 气象 因子 下 CAR 模型 构建 与 验证 

将 锡林郭勒 盟 1980 一 2019 年 的 气象 数据 集 按 
时 间 划 分 为 率 定期 与 验证 期 ,其 中 1980 一 2009 年 为 
率 定期 率 定 模型 ,2010 一 2019 年 为 验证 期 ,验证 模 
型 精度 。 按 照 PM-ET, 与 各 气象 因子 关联 系数 从 大 
到 小 的 顺序 建立 CAR 模型 ,模型 所 用 的 递 推 最 小 乘 
法 的 遗忘 因子 都 为 1.0, 模 型 的 显著 水 平 为 0.05 , 检 
验 结果 为 F=2.6896。 其 中 CAR, 为 以 关联 系数 从 大 
到 小 排列 的 气象 因子 数 的 建立 的 模型 ,n 为 气象 因 
子 数 。 通 过 比较 CAR-ET, 模 型 与 PM-ET, 模 型 的 拟 
合 关 系 ,分 析 .评价 所 建立 的 CAR-ETo 模 型 的 建 模 效 
果 ( 图 3、 图 4)。 

CAR-ETo 模 型 在 率 定 期 与 PM-ET。o 拟 合 较 好 ,在 
验证 期 稍 差 ,而 PM-ET, 与 HS-ET, 在 各 气象 站 点 预测 
值 均 低 于 PM-ET,, 拟 合 精度 较 差 。 模 型 在 率 定 期 内 
随 着 输入 气象 因子 的 数量 增加 建立 模型 的 尼 增 加 ， 
验证 期 内 模型 尼 随 着 输入 气象 因子 的 增加 呈现 先 
增加 后 减少 的 变化 趋势 ,表明 随 输 入 参数 增多 模型 
预测 值 更 接近 标准 值 。 模 型 的 NSE 在 模型 率 定期 
与 模型 验证 期 除 东 乌 珠 穆 沁 旗 外 都 在 0.9 以 上 , 且 


各 CAR-ETo 模 型 的 NSE 值 都 大 于 HS-ET, 的 NSE 值 ， 
RMSE 同 尼 也 有 着 相似 的 变化 规律 ,在 模型 率 定 期 
随 着 输入 气象 因子 的 增加 ,CAR-ET, 与 PM-ET。 计 算 
值 的 RMSE 减 小 ,验证 期 内 ,模型 的 RMSE 值 并 没有 
随 着 输入 的 气象 因子 增加 而 减 小 

通过 比较 建立 的 CAR 模 型 与 HS-ET,、PMT-ET。 
只 需要 气温 数据 简化 计算 方法 计算 精度 来 确定 模 
型 , 阿 巴 嘎 着 确定 为 CAR; 模 型 ,其 RR、NSE、RMSE 参 
数 均 优 于 HS-ET, 与 PMT-ET, 简 化 计算 方法 , 东 乌 珠 
穆 沁 旗 为 CAR: 模 型 , 除 尼 值 为 0.646 小 于 HS-ET 方 
法 0.681 外 ,其 他 参数 值 均 优 于 2 个 简化 计算 方法 ， 
多 伦 县 为 CAR: 模 型 , 除 尼 值 为 0.392 小 于 HS-ET 方 
法 0.782 外 ,其 他 参数 值 均 优 于 2 种 简化 计算 方法 ， 
二 连 浩特 为 CAR; 模 型 ,所 有 检验 参数 都 优 于 PMT- 
ET 与 HS-ET, 算 法 , 西 乌 珠 穆 沁 旗 为 CARs 模 型 , 除 
尼 值 为 0.614 小 于 HS-ET" 方 法 0.739 外 ,其 他 参数 值 
均 优 于 PMT-ET, 与 HS-ETo 算 法 ,锡林浩特 为 CAR; 模 
型 , 除 尺 值 为 0.767 小 于 HS-ET。o 与 PMT-ET, 计 算 方 
法 的 尼 值 0.826 .0.829 外 ,其 他 参数 全 都 优 于 PMT- 
ET 与 HS-ETu 算 法 。 通 过 上 述 比 较 得 出 ,处 于 草 多 
草原 区 的 东 乌 珠 穆 沁 旗 与 西 乌 珠 穆 沁 旗 仅 需 日 昭 
时 数 与 最 高 气温 数据 2 个 气象 因子 数据 就 可 以 较 高 
精度 预测 出 参考 作物 蒸 散 量 有 7v 的 值 ,处 于 典型 草 
原 区 的 阿 巴 嘎 旗 锡林浩特、 多 伦 县 需要 除 日 照 时 
数 .最 高 气温 数据 外 还 需要 平均 风速 数据 与 平均 相 
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a= PM-ETo == PMT-ET, =h HS-ETo i CAR, =0= CAR; 
(a) 阿 巴 嘎 旗 b) 东 乌 珠 穆 沁 旗 
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图 3 CAR 模 型 在 模型 率 定 期 和 验证 期 内 与 PM-ET 的 拟 合 效果 
Fig. 3 The fitting effect of CAR model with PM-ET during model rate and validation period 
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注 : 


C, 表 示 CAR 模型 ," 为 气象 因子 个 数 ,再 表示 HS-ETu,P 表 示 PMT-ET。 


图 4 各 模型 评价 参数 比较 


Fig. 4 Comparison of evaluation parameters of each model C, represents CAR model, n represents the number 


of meteorological factors, H represents HS-ET», P represents PMT-ET) 
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对 湿度 数据 才能 以 较 高 精度 预测 E76 的 值 , 处 在 荒 
漠 草 原 区 的 二 连 浩特 仅 需 日 照 时 数 与 最 高 气温 数 
据 2 个 气象 因子 数据 就 可 以 较 高 精度 预测 出 参考 作 
PAR BCE ET) WEL. 


3 讨论 


目前 对 参考 作物 蒸 散 量 的 预测 较为 通用 的 方 
法 为 PM-ET, 方 法 ,但 能 够 提供 ET, 计算 所 需要 的 全 
部 气象 因子 的 国家 气象 站 均 数 不 足 1 个 ,这 给 7 
精准 预报 普及 带 来 了 很 大 困难 。 在 国家 级 气象 站 
数据 支持 下 ,可 以 根据 PM-ET, 模 型 计算 出 87 标准 
值 ,利用 国家 级 气象 站 辐射 范围 以 外 的 区 域 级 气象 
站 (能 提供 87, 计 算 所 需 部 分 气象 因子 ,但 数目 远 多 
于 国家 级 气象 站 、 建 设 费 用 远 低 于 国家 级 气象 站 ) 
数据 ,构建 模型 较为 准确 的 对 ET 进行 预报 ,对 气候 
干 湿 状 况 分 析 、 区 域 水 资源 评价 与 农田 水 分 管理 与 
灌溉 决策 意义 重大 。 

由 于 锡林郭勒 主要 处 在 干旱 、 半 干旱 气候 区 ， 
有 着 风速 大 且 风 速 在 空间 与 时 间 的 强 变 异性 特殊 
的 气候 特征 ,导致 HS-ET, 与 PMT-ET, 在 不 考虑 风速 
情况 下 均 使 pT 被 低估 (图 3), 同 时 难以 得 到 合适 的 
风速 与 风速 变异 数学 规律 ,使 得 上 述 2 种 方法 对 ET 
预测 精度 提高 较为 困难 下 。 多 变量 时 间 序 列 CAR 
模型 是 在 只 有 一 个 变量 的 积分 自 回 归 滑 动 平 均 
(ARIMA ) 模 型 基础 上 改进 得 到 的 "3, 模型 具有 一 维 
时 间 序 列 与 回归 分 析 2 类 数理 统计 方法 的 优点 , 考 
JENA AFM Smo” ,符合 参考 作物 蒸 散 量 的 变化 
规律 。 本 文选 取 了 锡林郭勒 盟 的 国家 级 地 面 气象 
站 点 ,通过 关联 度 分 析 得 到 了 各 气象 因子 与 PM-ET。 
关联 系数 ,应 用 反 距离 插值 法 得 出 关联 系数 在 研究 
区 的 空间 分 布 “-, 按 与 多 年 平均 PM-ET。 的 关联 系 
数 从 大 到 小 的 顺序 选择 气象 因子 建立 多 变量 时 间 
序列 CAR-ET, 模 型 ,将 PM-ET, 设 为 标准 方法 ,再 与 
HS-ET, 法 进行 比较 ,评价 模型 。 研 究 表明 ,CAR-ET。 
模型 结合 特定 象 因子 数据 的 情况 下 预测 精度 较 HS- 
ET 与 PM-ET, 方 法 更 高 ,但 是 模型 均 要 求 输入 日 照 
时 数 等 与 PM-ET, 关 联 度 较 大 气象 参数 ,而 日 照 时 数 
又 较 难 获取 ,需要 当地 依照 实际 制定 数据 监测 
方案 。 

锡林郭勒 草原 面积 广阔 ,然而 国家 级 的 地 面 气 
象 站 辐射 的 范围 有 限 , 同 时 由 于 其 地 理 差 异 使 得 参 


ChinaXiv 合 作 期 刊 


6 期 冯 壮 壮 等 :基于 多 变量 时 间 序列 模型 的 锡林郭勒 草原 参考 作物 蒸 散 量 预测 1655 


考 作物 蒸 散 量 在 锡林郭勒 草原 不 同 区 域 受 不 同 的 
气象 因子 影响 ,使 锡林郭勒 难以 使 用 PM-ETo 法 精确 
预测 参考 作物 蒸 散 量 ,可 通过 利用 设置 在 旗 ( 县 、 
市 .区 ) 级 行政 区 域 的 地 方 气象 观测 站 和 牧 业 气象 
试验 站 的 气象 数据 以 及 有 针对 性 布设 气象 监测 设 
备 ,收集 数据 构建 CAR-ETv 模 型 来 精准 预测 ET, AT 
进一步 为 锡林郭勒 草原 生态 环境 恢复 过 程 及 其 科 
学 问题 提供 数据 支撑 。 


4 结 论 

(1) 在 研究 区 内 ,各 气象 因子 与 PM-ET, 的 关联 
系数 在 研究 区 呈现 由 荒漠 草原 区 一 典型 草原 区 
草 甸 草 原 区 过 渡 减 小 的 分 布 ,其 中 日 最 低 气 温 HE 
均 气 温 与 PM-ET 关 联 度 最 小 ,平均 日 照 时 数 与 PM- 
ET 关联 度 最 大 ,为 0.7293 。 

(2) CAR-ET 模 型 除了 尼 小 于 HS-ET, 与 PMT- 
ET 模型 外 ,由 于 HS-ET, 与 PMT-ET, 模 型 对 ETv 低 
估 , CAR- ET, 预测 精度 总 体 优 于 HS- ET, 45 PMT- 
ET,。 处 于 锡林郭勒 草 多 草原 区 .荒漠 草原 区 通过 日 
照 时 数 、 最 高 气温 二 因子 建 模 计算 得 到 的 87, 精 度 
最 高 ,锡林郭勒 典型 草原 区 通过 日 照 时 数 .最 高 气 
温 、 风 速 三 因子 据 建 模 计算 精度 最 高 。 

(3) 锡林郭勒 草原 受 其 特殊 的 地 理 条 件 与 特殊 
气候 影响 ,只 通过 输入 气温 参数 的 PMT-ETo 与 HS- 
ET 方法 精度 低 于 考虑 日 照 . 风 速 影响 的 CAR-ETo 模 
型 ,可 结合 当地 地 理 与 气候 特征 因地制宜 制定 气象 
监测 方案 建立 CAR-ET, 模 型 进行 57 预测。 
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Prediction of reference crop evapotranspiration in Xilinguole grassland based 


on multivariate time series model 


FENG Zhuangzhuang, SHI Haibin, MIAO Qingfeng, LI Jiannan, 
SUN Wei, DAI Liping 
(College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, 


Hohhot 010018, Inner Mongolia, China) 


Abstract: This study was performed to explore the influence of various meteorological factors on reference crop 
evapotranspiration and its prediction in Xilinguole grassland. The correlation between the meteorological data of 
six national ground meteorological stations in Xilinguole grassland of Inner Mongolia and the calculated PM-ET) 
was analyzed. A multivariate time series controlled autoregressive (CAR) model was used as a basis for 
establishing a CAR-ET» model in accordance with the order of correlation coefficient. Results showed that the 
correlation coefficient between meteorological factors and PM-ET decreased from a desert steppe area to a typical 
steppe area and a meadow steppe area. The correlation of daily minimum temperature, average temperature, and 
PM-ET。 was the smallest, whereas the correlation between average sunshine hours and PM-ET) was the largest, 
which was 0.7293. In meadow and desert steppes, the accuracy of CAR-ET) increased when sunshine and temperature 
were used as inputs. Wind speed was also required as the third factor of the model in a typical grassland area. The 
verification of model prediction accuracy revealed that the prediction accuracy of the CAR- ET» model was 
generally higher than that of HS- ET, and PMT- ET) models. Combined with correlation analysis, the 
meteorological factors with a significant impact on local reference crop evapotranspiration could be obtained. 
Thus, this study could provide a theoretical basis for designing a meteorological monitoring layout, determining 
reference crop evapotranspiration, and conducting ecological restoration in Xilinguole grassland. 

Keywords: reference crop evapotranspiration; grey correlation analysis; multivariable time series model; Xilin- 


gol grassland 


